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Chiral quark models with non-local covariant separable interactions at finite temperature 
and chemical potential are investigated. We develop a formalism in which the different quark 

(Zj ■ properties are evaluated taking into account the analytic structure of the quark propagator. 

^ ' In this framework we study the chiral restoration phase transition for several definite non- 

local regulators, including that arising within the instanton liquid picture. We find that in 
Cn ■ all cases the chiral transition is of first order for low values of T, turinng into a smooth 

crossover at a certain "end point" . Using model parameters which lead to the physical pion 
mass and decay constant, we find for the position of this "end point" the values (Tb, /i_b) ^ 
(60 — 70, 180 — 210) MeV. We also discuss the special relevance of the first poles of the quark 
propagator. 
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I. INTRODUCTION 
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I I The understanding of the behaviour of strongly interacting matter at finite temperature and/or density 

Jr~^' is of fundamental interest and has important applications in cosmology, in the astrophysics of neutron stars 

(-H ' and in the physics of relativistic heavy ion collisions. From the theory of the quark-gluon interactions. 

Quantum Chromo Dynamics (QCD), we beheve that at zero temperature and density chiral symmetry is 
spontaneously broken and quarks are confined within hadrons. However, since QCD is asymptotically free, 
k> , when either the temperature T or the chemical potential /i are high the effective coupling becomes small. 

^ • Thus, one expects that at a certain point the system undergoes a phase transition (or crossover) to a new 

^ I phase in which color is screened rather than confined and chiral symmetry is restored. Unfortunately, so far, 

it has not been possible to obtain detailed information about the corresponding T — fi phase diagram directly 
from QCD. In fact, lattice simulations, which work well for zero density and finite temperature, have serious 
difhculties in dealing with the complex fermion determinant that arises at finite chemical potential [y. In 
this situation, different models have been used to study this sort of problems. Among them the Nambu- 
Jona-Lasinio model ^ is one of the most popular. In this model the quark fields interact via local four 
point vertices which are subject to chiral symmetry. If such interaction is strong enough chiral symmetry is 
spontaneously broken and pseudoscalar Goldstone bosons appear. It has been shown by many authors that 
when the temperature and/or density increase, the chiral symmetry is restored pi . As an improvement on the 
local NJL model, some covariant non-local extensions have been studied in the last few years |j]. Nonlocality 
arises naturally in the context of several of the most successful approaches to low-energy quark dynamics 
as, for example, the instanton liquid model [pj and the Schwinger-Dyson resummation techniques 0. It 
has been also argued that non-local covariant extensions of the NJL model have several advantages over the 
local scheme. Indeed, non-local interactions regularize the model in such a way that anomalies are preserved 
||7| and charges properly quantized, the effective interaction is finite to all orders in the loop expansion and 
therefore there is not need to introduce extra cut-offs, soft regulators such as Gaussian functions lead to small 
NLO corrections |8|, etc. In addition, it has been shown Q that a proper choice of the non-local regulator 
and the model parameters can lead to some form of quark confinement, in the sense that the effective quark 
propagator has no poles at real energies. In a previous letter llQ] we have studied the phase diagram of a 



tPellow of CONICET, Argentina. 



non-local model with a Gaussian regulator. The aim of the present work is to present details of such analysis 
as well as to extend it to more general cases, such as the Lorentzian regulator and the regulator arising 
within the instanton liquid model. 

This article is organized as follows. In Sec. II we introduce the model and the formalism needed to 
extend it to finite temperature and density. In Sec. Ill we discuss how to perform the summation over the 
Matsubara modes for the expressions obtained previously. Our results for some specific non-local regulators 
are presented in Sec. IV. Finally, in Sec. V we give our conclusions. 

II. FORMULATION 

Our starting point is the partition function at zero T and /^, 

Zo ^ J V^V^, e-^^ , (1) 

where Se stands for the Euclidean action 

Se — I d'^x 



G 

ll^{x) {-ip + nic) Ipix) - —ja{x)jaix) 

Here rric is the current quark mass, and the euclidean operator ^ is defined as 



(2) 



^ = 74|^+7-V (3) 

with 74 = z7o, r = it. The current ja{x) is given by 

ja{x) = / d*y d*z r{y - x) r{x - z) ■i}){y)Tai}{z) , (4) 



where Ta = (1, ^75^) and r{x — y) is a non-local regulator function. The regulator is local in momentum 
space, namely 

^{^-^)- j ^.^-'^'^'^"rip). (5) 

In fact, Lorentz invariance implies that r{p) can only be a function of p^, hence we will use for the Fourier 
transform of the regulator the form r(p^) from now on. 

To proceed it is convenient to deal with bosonic degrees of freedom. Let us perform a standard bosonization 
of the theory, introducing the sigma and pion meson fields 

Ma{x)^{a{x),Ti{x)) . (6) 

Following the usual steps we obtain a partition function equivalent to that in (fil): 



VaVi: delA{Ma) exp ^— j -^ M^p) 



(7) 



where the operator A in momentum space reads 

A{Ma) = i-^ + m,){27r)U^''Hp-p') + r{p')Maip-p')r{p'^)Ta . (8) 

At this stage, we perform the mean field approximation by expanding the meson fields around their trans- 
lational invariant vacuum expectation values 

(j(x) = CT+ 5a{x) (9) 

TTi (x) = vfj -f- S7ri{x) (10) 

and neglecting the fiuctuations 6a (x) and S'!Ti{x). The mean values of the pion fields vanish due to symmetry 
reasons. Within this approximation the determinant in (H) is formally given by 



/^4 
tr log [-i> + mc + a r^ [p^) ] 
(ZTTJ 



(11) 



where tr stands for the trace over the Dirac, flavor and color indices, and V^^' is the four-dimensional volume 
of the functional integral. 

Now the corresponding partition function in the grand canonical ensemble for finite temperature T and 
chemical potential ji can be obtained through the replacement in the integrals in ([7|) and (|ll| ) 

[^F{p,,p) ^ ^F{p,,p)^T^-^ f^Fiu;^-^f,,p)^ (12) 

where uin are the Matsubara frequencies corresponding to fermionic modes, u)n — (2n + 1)'kT . In performing 
this replacement we have assumed that the regulator acquires an explicit /i dependence, as it is the case e.g. 
in the instanton liquid model [ [ll]Jl^ . In the same way we replace the volume V''^^ by V/T, V being the 
three-dimensional volume in coordinate space. The mean field partition function reads then 

Zmf{T, /i) = exp [-Tmf{T, ^)] , (13) 

where the mean field action is given by 

|rMf(r,M) = -47vJ^iog[p2 + s2(p2)]j +g (14) 

with the quark selfenergy S(p^) defined as 

S](p2)=m, + ar2(p2). (15) 

Finally, the grand canonical thermodynamic potential is given by 

nMF{T,^x) = ~y log ZMF{T,^i) = ^TMFiT,fl). (16) 



It can be shown that in general this quantity turns out to be divergent. However, it can be regularized 

^MF ~ ^'■MF — ^'■MF- 



by subtracting the corresponding value at zero T and /i, f^»Vp = ^mf — ^mf- ^^ow the minima of the 



thermodynamic potential are obtained from the solutions of 

dQ.MF 



which leads to the following gap equation for a: 

a-m,G^ p2 + s2(p2) ■ (18) 

Given the thermodynamic potential the expressions for all other relevant quantities can be easily derived. 
For each flavor the quark vacuum expectation value is given by 

while the quark density turns out to be 

In obtaining these last two equations it should be noticed that for each quark flavor only half of the kinetic 
term in Eq.dlJ) contributes to the derivatives. 



As well known, if the quark selfenergy is momentum independent as in the conventional NJL model, 
the sums over the Matsubara frequencies indicated in the previous equations can be easily carried out and 
one obtains rather simple expressions in terms of the usual occupation numbers. This also holds when the 
selfenergy only depends on the spatial components of the momentum. However, for the covariant regulators 
we are considering here the situation is more complicated. The main difficulty is that the analytic structure 
of the quark propagator S{q) = l/(— ^ + S(q^)) in the complex plane can be much richer in this case: there 
might be a rather complicate structure of poles and cuts. 

In Minkowski space and for q — Q the positions of the poles are given by the solutions of 

qo^-[mc + ar^i-qo^)f=0. (21) 

Let us denote the real and imaginary parts of these solutions with rp and ip respectively. For nonzero spatial 
momentum (fit is easy to see that the poles are located at Zp = Rp + Up, where 

Rp = ±ep , Ip = ±^^ (22) 



e 
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with 



^p-ip + 9^ + V(rp"ip + 9^)^+4 r-2 ij 



'P^\ 2 ■ ^^^^ 

On the other hand, the possible cuts of the propagator will be given in general by the cuts of the regulator 
as a function of q. As in the case of the poles, it is convenient to determine first their position in the qo 
plane for q = 0. Then it is simple to find the position in the general case applying the above equations to 
all the points along the cuts. 

Before going into the explicit evaluation of the Matsubara sums it is useful to analyze the pole and cut 
configurations in some relevant situations. One important case is that in which the regulator is such that 
the Minkowski quark propagator has an arbitrary set of poles but no cuts. Within this class of regulators it 
might exist a situation in which there are no poles along the real axis. As already mentioned, this situation 
might be interpreted as a realization of confinement [p|. In that case only quartets of poles located at 
ap = Tp ± i ip, ap = —rp ± i ip appear. On the other hand, if purely real poles do exist they will show 
up as doublets ap — ± rp. It is clear that the number and position of the poles depend on the details of 
the regulator. For example, if we assume it to be a step function as in the standard NJL model, only two 
purely real poles at ± M appear, with M being the dynamical quark mass. For a Gaussian interaction, three 
different situations might occur. These can be classified according to the value of S(0) at zero T and /i, 
which we denote by S(0). For values of S(0) below a certain critical value l](0)crit, two pairs of purely real 
simple poles and an infinite set of quartets of complex simple poles appear. At S(0) = S(0)crit, the two pairs 
of purely imaginary simple poles turn into a doublet of double poles with ip = 0, while for S(0) > T,{0)crit 
only an infinite set of quartets of complex simple poles is obtained. In the case of the Lorentzian interactions, 
there is also a critical value above which purely real poles at low qg cease to exist. However, for this family 
of regulators the total number of poles is always finite. As we see, several physically interesting regulators 
indeed belong to the "no cuts" group. Unfortunately some other important ones, like that arising within 
the instanton liquid model, do lead to propagators which present poles and cuts. Thus, it is necessary to 
consider this more general situation. 

III. EVALUATION OF MATSUBARA SUMS 

To proceed with the evaluation of the sums over the Matsubara frequencies we assume that the regulator is 
such that the quark propagator at (f = has the analytic structure shown in Fig. la. Whereas the proposed 
pole distribution is completely general, we consider for simplicity the particular case in which there is a 
single cut lying along the real axis in the regions [— oo, —rd and [re, cxd]. Let us consider the sum appearing 
in the gap equation Eq. (p^ , 



^- / (03 ^ E F[i^n-^,f + q')]=J^I{q') , (24) 



where 



F{q' 



g2 + S2(g2) 



(25) 



In order to evaluate the sum in (24) it is convenient to fix (f and define J^{z) = F\{—iz — i/i)^ -^ q"^)], so that 
the argument of the sum is given by JF(«a;„). It is easy to see that the structure of cuts and poles of J^{z) is 
similar to that shown in Fig. la, conveniently translated to the z plane (see Fig. lb). This translation can 
be performed by means of Eqs.(22), with an additional shift Rp —t Rp — /^. We introduce now the auxiliary 
function 



/(z) = [1 + exp(z/r)]- 



(26) 



which has a set of poles at z„ = iujn with — T as the corresponding residues. Provided that J^{z) has no 
poles on the imaginary axis, the sum we want to evaluate can be regarded as the sum over the z„ poles of 
the residues of the function — JF(z)/(z). We can use Cauchy's theorem to get 



I{<1 



T J2 n^ 



n— — oo 
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■zoo+(5 



dz T{z) f{z) 



— ZCXD + (5 



2tt 
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dz T{z) fiz) 



(27) 



— ioo — (5 



The integrals in Eq. ( p7| ) can be calculated with the help of complex integrals along adequate paths in the z 
plane. Let us consider the sum 

dzT{z)f{z) + I dzT{z)J{z) + / dzT{z)f{z) , (28) 

where the contours Ci, C2 and C3 are shown in Fig. lb. In the limit i? — > 00, one can make use of Cauchy's 
theorem to obtain 

I{q^)= Y. ^cs[T{z):Zp]f{zp) 
1 



dx AJ^{x)f{x) 



27r 



where Zp = Rp + Up are the poles of ^{z), and we have defined 

^{x + ie) — J-{x — ie) 



dz T{z) f{z) 



AT: 



2i 



(29) 



(30) 



which is different from zero along the cut. Notice that for large R the integrals over the round contours are 
expected to vanish in view of the exponential fall of f{z). 
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FIG. 1. (a) Poles and cuts of the quark propagator in the go plane for q — 0. (h) Pole and cut structure of the 
oropagator in the z plane and useful integration paths to evaluate the Matsubara sums. 



Now the last integral in the right hand side of (E3) can be rewritten by considering the function 

loiq^) = -7r dtz T{z) = — / d<Z4 F(g2) , (31) 

^^ J-ioo-t, ^^ J-oo 

which is nothing but I{q'^) in the limit of vanishing T and /i. Notice that 

/zoo — /i - pioo — fi 

dzT{z)J{z) = h[q^) + - dzT{z)f{-z). (32) 

-ioo — fi '^'^ J -ioo-fi 

The last integral in this expression can be evaluated by closing a contour on the left side of the complex 
plane and using once again Cauchy's theorem. Finally, assuming that the regulator is an analytic function 
of q^ , we can use the properties 

AJ^{—x — ^) ~ —lS.J-{x — /i) == — Im J-{x — fi + ie) 

Res [J^{z); Zp] — —Res [J-{z); —Zp — 2fi] , (33) 

to end up with our main result 

^^/(^^('Z') + Sy^Re Y. 7pRes[.F(z);Zp] [n+(zp + /i)+n_(zp + /i)] 



Rp > -/J 
/„ > 



1 f (fq 



'^ J (27r)3 y^7fTi= 



dx Im !F{x ~ II + ie) [nj^{x) + n_(x) ] , (34) 



where the first term corresponds to the value of 5* at zero temperature an chemical potential. The coefficient 
7p is defined as jp = 1/2 for Ip — and 7p = 1 otherwise, and n±{z) stand for the occupation number 
functions 

1 + exp [{z T w/rj 

It is clear that the steps leading to Eq.(E3) can be also followed to evaluate the sums appearing in e.g. the 
quark vacuum expectation value and the quark density, Eqs.dl9) and (EQ) respectively, once the function 
F{q^) is properly redefined. In the case of the regularized thermodynamic potential ^j^Jp the situation is 
somewhat more complicated, since the argument of the sum includes a logarithm that introduces cuts in the 
z plane outside the real axis. Nevertheless a similar relation can be derived to calculate the corresponding 
Matsubara sum. 

From Eq.(p4), it is now easy to obtain the well-known results for the standard NJL model with a three- 
momentum cut-off A ||3|. In fact, at the level of mean field approximation such model can be viewed as a 
particular case of the non-local schemes analyzed above, where the regulator is taken as a theta function in 
three-momentum space, r{q^) = 9{q^ — A^) (i.e. it is not covariant). Notice that for this regulator there is 
only one pole at zq = y/q"^ + M^ — u, M — mc + a, with Res[J^{z); zq] — — M/2(zo -I- ^). From the relation 
S = (t/{8GNc) (see Eqs.(p4h and (ttSh) one immediately arrives at the standard NJL gap equation 

2GN P^ M 

M^m, + ^ / dpp^ —[1^ n+ (Ep) - n_ (Ep)] , (36) 

'"' Jo ^p 



where Ep = yJ^^^^TwP- . 

IV. RESULTS FOR DEFINITE NON-LOCAL REGULATORS 

Having introduced the formalism needed to deal with models with non-local separable interactions at finite 
temperature and chemical potential, we turn now to show the results of numerical calculations carried out 
for some particular regulators. We concentrate here in three cases: the Gaussian regulator, the Lorentzian 
regulator and the regulator arising within the instanton liquid model. 



A. Gaussian regulator 

Let us analyze a regulator function of the form 

r(g2) = cxp(-gV2A2) , 



(37) 



where A is a parameter of the model, besides the coupling constant G and the current quark mass rric- We 
consider here two sets of values for the parameters. Set I corresponds to G = 50 GeV~^, rric = 10.5 MeV 
and A = 627 MeV, while for Set II the respective values are G = 30 GeV"^, rric = 7.7 MeV 
and A — 760 MeV. Both sets of parameters lead to the physical values of the pion mass and decay 
constant. At zero temperature and chemical potential, the calculated values of the chiral quark condensate 
are —(200 MeV)"^ for Set I and —(220 MeV)"^ for Set II. These values are similar in size to those determined 
from lattice gauge theory or QCD sum rules. The corresponding results for the selfenergy at T = /i = and 
zero momentum are f](0) = 350 MeV for Set I and i;(0) = 300 MeV for Set II. 

As stated in Sect. II, in the case of a Gaussian regulator the quark propagator has no purely real poles if 
E(0) is above a critical value S(0)crit, whereas for S(0) < S](0)crit two pairs of real poles exist. From the 
explicit expression for the critical selfenergy. 



S(0)c 



l{V^c 



2A2 - TOc ] exp - ( ymfT2A2 



/4A^ 



(38) 



it can be easily checked that Set I and Set II correspond to the first and second situations, respectively. Thus, 
following Ref. Q, Set I might be interpreted as a confining one, since quarks cannot materialize on-shell in 
Minkowski space. 
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FIG. 2. Behaviour of the selfenergy (upper graphs) and the density (lower graphs) with the chemical potential for 
three representative values of the temperature. Full lines corresponds to T — 0, dashed lines to T — 50 MeV and 
dotted lines to T = 100 MeV. The left panels display the results for Set I and the right panels those for Set II. 



The behaviour of the zero- momentum selfenergy S(0) as a function of the temperature and chemical 
potential can be obtained by solving the gap equation for finite T and /i, Eq.(p8|), where the Matsubara 



sum can be evaluated using Eq.(p3). In this case the regulator has no cuts in the complex plane, hence 
the last term of ( p4|) v anishes. A similar procedure can be followed to obtain the quark density p for finite 
T and /x from Eq. (|20|) . Our numerical results are shown in Fig. 2, where we plot S(0) vs. /x and p vs. pi 
for some fixed values of the temperature. The values of p are given with respect to nuclear matter density, 
po — 1-3 X 10^ MeV'^. The left and right panels in the figure correspond to the results for Set I and Set II, 
respectively. For both sets of parameters we observe the existence of some kind of phase transition at (or 
around) a given value of the chemical potential which depends on the temperature. The only qualitative 
differences appear in the behaviour of S(0) and p as functions of p for low temperatures. For example, 
for T — hO MeV we observe that in the case of the confining set (Set I) the selfenergy rises with /i (for 
p below the critical chemical potential) while the opposite occurs in the case of Set II. As one increases 
the temperature the behaviour of Set I turns into that of Set II and, eventually, in both cases we obtain a 
continuous and always decreasing curve. An equivalent behaviour is observed in the case of the density. To 
obtain our results, we have included in the sums over the poles of the quark propagator only the first few 
poles close to the origin in the go plane. We have checked, however, that for the range of values of T and 
p covered in our calculations, the convergence is so fast that already the first pole gives almost 100% of the 
full result. Thus, the behaviour of the relevant physical quantities up to (and somewhat above) the phase 
transition is basically dominated by the first set of poles (doublet or quartets) of the propagator. 

We observe that at T = there is a first order phase transition for both the confining and the non-confining 
sets of parameters. As the temperature increases, the value of the chemical potential at which the transition 
shows up decreases. Finally, above a certain value of the temperature the first order phase transition does 
not longer exist and, instead, there is a smooth crossover. This phenomenon is clearly shown in the left 
panel of Fig. 3, where we display the critical temperature at which the phase transition occurs as a function 
of the chemical potential. The point at which the first order phase transition ceases to exist is usually called 
"end point" . 
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FIG. 3. Critical temperatures as a function of the chemical potential. The left panel corresponds to the case of finite 
quark masses and the right panel to the chiral limit. In the right panel, the dotted lines correspond to the region of 
the second order phase transition, while the full lines in both panels indicate the region where a first order transition 
occurs. The dashed line in the right panel indicates the critical line for deconfinement corresponding to Set I. For 
chemical potentials somewhat larger than pp this line coincides with that of the chiral restoration. 



In the chiral limit, the "end point" is expected to turn into a so-called "tricritical point", where the second 
order phase transition expected to happen in QCD with two massless quarks becomes a first order one. 
Indeed, this is what happens within the present model for rric = 0, as it is shown in the right panel of Fig. 3. 
Some predictions about both the position of the "tricritical point" and its possible experimental signatures 
exist in the literature |13|lJ]. In our case this point is located at {Tp,pp) =(70 MeV, 130 MeV) for Set I 
and (70 MeV, 140 MeV) for Set II, while the "end points" are placed at {Te,pe) =(70 MeV, 180 MeV) and 
(55 MeV, 210 MeV), respectively. 

In general, for temperatures below the end point there is a range of values of p where one finds two solutions 
of the gap equation, i.e. there are two values of the selfenergy where the thermodynamic potential shows 
a local minimum. The situation is illustrated in Fig. 4, where we plot the (renormalized) thermodynamic 
potential vs. the selfenergy for different values of the temperature and some fixed values of the chemical 
potential around the phase transition point. The curves in the figure correspond to the parameters of Set II. 



In the left panel we consider the case of T = 0, where the phase transition occurs at Hc — 270 MeV, and 
we plot the values of fi^'"^ for values of fi which are 10 MeV above and below Hc- The equilibrium point 
is clearly that corresponding to the absolute minimum in each case. In the central panel we take T — 50 
MeV, close to the end point temperature of ^ 55 MeV, hence the two minima approach each other and the 
thermodynamic potential between them is approximately constant. For the right panel the temperature is 
100 MeV, well above the end point, therefore one finds only one minimum and the phase transition proceeds 
through a smooth crossover. 
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FIG. 4. Renormalized thermodynamic potential for Set II as a function of the selfenergy, for different values of the 
chemical potential fi. Left, central and right panels correspond to T = 0, T = 50 MeV and T — 100 MeV respectively. 

The behaviour of the critical temperatures with the quark density is shown in Fig. 5, where we plot the T—p 
phase diagrams for both the confining and the non-confining sets. The curves in the left panel correspond to 
the case of finite quark masses, and the fat dots indicate the end points for each set of parameters. The area 
below each curve is a region where both phases are allowed. This "mixed phase" zone can be interpreted |13t| 
as a region where droplets containing light quarks of mass rUc coexist with a gas of constituent, massive 
quarks. The corresponding curves in the chiral case {nic = 0) for both Sets I and II are shown in the right 
panel of the figure. Here the dotted lines correspond to the second order transitions, and the fat dots indicate 
the tricritical points. 
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FIG. 5. Critical temperatures as functions of the relative density for Sets I and II. The left panel corresponds to 
the case of finite quark masses and the right panel to the chiral limit. The fat dots indicate the end points (left panel) 
and tricritical points (right panel), and the dotted lines in the right panel indicate the second order phase transitions. 



It is interesting to discuss in detail the situation concerning the confining set. In this case we can find, 
for each temperature, the chemical potential /i^ at which confinement is lost. Following the proposal of 
Ref. 0, this corresponds to the point at which the selfenergy at zero momentum reaches the value E(0)crit- 
Using the values of rUc and A corresponding to Set I we get S(0)crjt — 267 MeV. We see that for low 
temperatures the value of fid coincides with the chemical potential corresponding to the chiral restoration, 
while for temperatures close enough to that of the "end point" it starts to lie slightly below the phase 



transition line. Above the end point temperature Te it is difRcult to make an accurate comparison since, 
for finite quark masses, the chiral restoration proceeds through a smooth crossover. However, we can still 
study the situation in the chiral limit. We find that in the region where the chiral phase transition is of 
second order deconfinement always occurs, for fixed T, at a lower value of /i than the chiral restoration. The 
corresponding critical line is indicated by a dashed line in the right panel of Fig. 3. In any case, as we can 
see in this figure, the departure of the line of chiral restoration from that of deconfinement is in general not 
too large. This indicates that within the present model both transitions tend to happen at, approximately, 
the same point. 



B. Lorentzian regulator 



Let us now consider the case of a Lorentzian regulator. 



r(g2) = 



1 



1 + (g2/A2)" ' 



(39) 



where n is a positive integer. For this family of regulators, it can be seen that the quark propagator always 
has a finite number of poles in the go plane. If n is even, there are 2n quartets of complex poles, plus a 
doublet of real poles. For odd values of n, one has a doublet (quartet) of real (complex) poles for S(0) below 
(above) a certain critical value l](0)crit, plus a set of 2?! — 1 quartets of complex poles, and a doublet of real 
poles. This last doublet occurs at Qq > A^. 

In order to analyze the characteristics of the chiral phase transition for this shape of the regulator, we 
have taken the simplest case n = 1, choosing as input parameter the quark selfenergy at zero temperature 
and chemical potential, S(0). As in the case of the Gaussian regulator, the parameters G, rric and A have 
been fixed so as to reproduce correctly the physical values of the pion mass and decay constant. For a 
typical selfenergy S(0) = 300_MeV we find the reasonable values G = 28.38 GeV"^ m^ = 4.57 MeV and 
A = 940 MeV. The value of l](0)crjt is in this case approximately 270 MeV, hence the fermion propagator 
has two quartets and one doublet of poles. This is represented in Fig. 6. 
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FIG. 6. Analytic structure of the quark propagator in the go plane for an n = 1 Lorentzian regulator. 



Once again, we can make use of both the gap equation and our result Eq. (pm to evaluate the quark self- 
energy S(0) for given values of the temperature and chemical potential. The corresponding phase transition 
curves are plotted in the left panel of Fig. 7, where we have considered fixed temperatures of 0, 30, 50 and 
100 MeV. As in the previous cases, we observe that for low values of T the chiral phase transition is of 
first order, while for values above a certain Te the transition proceeds through a smooth crossover. In the 
right panel of the figure we plot for the same values of the temperature the relative quark density p/po as a 
function of the chemical potential. 
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FIG. 7. Behaviour of the quark self energy (left) and the relative quark density (right) with the chemical potential in 
a non-local NJL model with a Lorentzian regulator, for some representative values of T . Full, dashed, dashed-dotted 
and dotted lines correspond to T = 0, 30, 50 and 100 MeV respectively. 

It is interesting to point out that for temperatures quite below Te (as e.g. 30 MeV for the present input 
parameters) and n < iic the selfenergy rises for increasing fi. Although the growth is too slight to be clearly 
perceived in the figure, it becomes evident at a magnified scale. Then, as T becomes closer to Te, this 
behaviour turns into a decreasing one (see the T — 50 MeV curve in the left panel of Fig. 7), ending up 
in the smooth crossover at T > Te- As mentioned in Sec. IV. a for the Gaussian interaction, the growing 
behaviour at low T is characteristic for the situation in which the Minkowski quark propagator has no real 
poles. We have stated above that for the Lorentzian regulator there is always a real pole close to gp ~ ^^• 
Thus, we can argue that the rise of T,(0) is mainly related to the absence of real poles at low ^q, which are in 
fact the poles to be associated with deconfined quarks. Moreover, this can be understood as a new evidence 
that the basic features of the T — /i behaviour of the quark properties are greatly determined by the first 
pole of the quark propagator. 

The values of the critical temperature as a function of the chemical potential and the quark density are 
displayed in the left and right panels of Fig. 8, respectively. Once again, the T — p phase diagram shows a 
mixed phase region. The predicted position of the end point is in this case {Te,He) =(58 MeV, 195 MeV). 
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FIG. 8. Critical temperature for the Lorentzian regulator as a function of the chemical potential (left) and the 
relative quark density (right). 



C. Instanton liquid model regulator 



The instanton liquid picture of the QCD vacuum predicts a separable non-local interaction with a regulator 
given by 



r{q^) = -x^ [Io{x)Ko{x) - h{x)Ki{x)] 



(40) 



x = Jq^pl2 



where / and K arc the modified Bessel functions and p stands for the average instanton size. Here we take 
the standard value p — 1/3 fm, and fix the couphng constant to G = 36.7 GeV~^ so as to reproduce the 
typical value for the instanton density n~ 1 fm^'' |^. Using these values, together with rric = 4.9 MeV, it is 
possible to reproduce the empirical values of the pion mass and decay constant. The resulting value of the 
chiral quark condensate at zero temperature and chemical potential is —(256 MeV)'^. 

As in the previous cases, we can solve the gap equation to obtain the phase transition curve for the quark 
selfenergy S(0). Here the situation is more complicated than for the Gaussian and Lorentzian cases, since 
the quark propagator has now a cut in the complex plane. As discussed in Sect. H, it is useful to look at the 
regulator function ( |40[) in Minkowski space for zero three-momentum, r{—q^). If we consider the analytic 
extension in the complex go plane, we see that the regulator shows a cut along the real axis, namely 



;^^21 



r{-{x + ief) -r{-{x-ief) = 2 i tt [2xJo{x)Ji{x) - Ji{x) 



(41) 



with X = px/2. It is easy to see that this cut is also present for the full quark propagator function. As 



stated in Sect. Ill, this translates into a cut on the real axis in the z plane if we set q 



-{z+^y 



■Q 



In 



addition, it can be seen that the quark propagator has a finite number of poles. The latter come in quartets, 
and their number increases with the value of the dimensionless parameter A = pS(0). One can find a set 
of critical values Xi, such that there is only one quartet of complex poles for < A < Aq, two quartets for 
Aq < A < Ai, etc. The analytic structure of the propagator for the case of two quartets is illustrated in 
Fig. 9. The critical values A^ only depend on the scaled current quark mass, prric. Taking p and nic as above 
we get for the first critical values Ao — 0.81, Ai ~ 1.09, A2 — 1.48. With the chosen value G — 36.7 GeV~^ 
the gap equation for T = p = leads to S(0) — 359 MeV, thus in our case A ~ 0.6 and we only have to 
consider a single quartet of complex poles. This means that the sum in Eq.(p4|) contains only one term. 
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FIG. 9. Analytic structure of the quark propagator for the instanton liquid model. 
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FIG. 10. Behaviour of the quark selfenergy (left) and the relative quark density (right) with the chemical potential 
for the instanton liquid model. Full, dashed, dashed-dotted and dotted lines correspond to T = 0, 50, 70 and 100 MeV 
respectively. 



The behaviour of the quark selfenergy and the quark density as functions of the chemical potential are 
shown in Fig. 10 for some representative values of the temperature. Once again, we find that the transition 
is of first order for low values of T, turning into a smooth crossover above a certain end point. As in the 
confining Gaussian and Lorentzian cases, we find that for finite temperatures quite below Te the selfenergy 
rises with /i for fi < ^c- Thus, this feature does not seem to be dependent on the existence of the cut, though 
we have checked that its contribution to the corresponding gap equation is not negligible. Finally, as in the 
previous cases, we present the curves corresponding to the phase transition temperatures as functions of /i 
and p/ po- These are shown in Fig. 11. For the instanton liquid model the end point is found to be located 
at {Te,pe) =(65 MeV, 180 MeV). 
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FIG. 11. Critical temperature for the instanton liquid model as a function of the chemical potential (left) and the 
relative quark density (right). 



V. SUMMARY AND CONCLUSIONS 



In this work we have studied the behaviour of chiral quark models with non-local covariant separable 
interactions at finite temperature and chemical potential. After introducing the general formalism we have 
derived some expressions that allow to compute the relevant quantities in terms of the analytic structure 
of the quark propagator. Although, in principle, the sum over the Matsubara frequencies could have been 
directly performed in a numerical way, we believe that our method allows for a better understanding of some 
features of the chiral restoration transition (as e.g. the importance of the first pole of the quark propagator) 
and stresses the connection with the usual procedure followed in the standard NJL model. In our numerical 
calculations we have considered three types of regulators: the Gaussian regulator, the Lorentzian regulator 
and the instanton liquid model regulator. In all these cases we have set the model parameters so as to 
reproduce the empirical values of the pion mass and decay constant and to get a chiral quark condensate in 
reasonable agreement to that determined from lattice gauge theory or QCD sum rules. We find that in all 
cases the phase diagram is quite similar. In particular, we obtain that for two light flavors the transition 
is of first order at low values of the temperature and becomes a smooth crossover at a certain "end point" . 
Our predictions for the position of this point are very similar for all the regulators and slightly smaller than 
the values in Refs. ||l|,|j, Te ~ 100 MeV and pE ~ 200 - 230 MeV. In this sense, we should remark that 
the model under consideration predicts a critical temperature at /.t = of about 100 MeV JT^ ] , somewhat 
below the values obtained in modern lattice simulations which suggest Tc ~ 140 — 190 MeV [Qpin any case, 
our calculations seem to indicate that pE might be smaller than previously expected even in the absence of 
strangeness degrees of freedom. 

An interesting feature of the behaviour of the quark selfenergy as a function of the chemical potential for 
temperatures well below Te is that it raises with p (for p < pc) when the first pole in the quark propagator 
is complex. A similar behaviour is found for the chiral condensate [|l^. This is the opposite to what 
happens when the first pole is real (in Minkowski space) as in the standard NJL model H or models where 
the nonlocality arises only in the space components |16[| . Lattice calculations, even at very low chemical 
potential as e.g. those reported in Ref. fl^, might be helpful to understand what happens in QCD. 

Several extensions of this work are of great interest. For example, it would be very important to investigate 
the impact of the introduction of strangeness degrees of freedom and flavor mixing on the main features of 



the chiral phase transition. In fact, some work on an SU{3) extension of the present type of models at finite 
/i, but T = 0, has been recently reported ||l^. Another topic that has recently attracted a lot of attention, 
and is suitable to be investigated in the present model, is the competition between chiral symmetry breaking 
and color superconductivity at large chemical potential. We expect to report on these issues in future 
publications. 
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